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We propose a realization of perfect valley filters based on the chiral domain-wall channels between 
a quantum anomalous Hall insulator and a quantum valley Hall insulator. Uniquely, all these 
channels reside in the same valley and propagate unidirectionally, 100% valley-polarizing passing-by 
carriers without backscattering. The valley index, the chirality, and the number of the channels are 
protected by topological charges, controllable by external fields, and detectable by circular dichroism. 


PACS numbers: 73.43.Cd, 71.70.Ej, 73.22.-f, 73.63.-b . 

Many crystalline materials have two or more degener¬ 
ate energy extrema in momentum space, which endows 
the low-energy carriers with a valley degree of freedom. 
Following the extensive studies in spintronics, there has 
been a growing interest in exploring the possibility of 
manipulating valley degrees for information processing, 
leading to a concept of valleytronics. 1-8 For practical val- 
leytronics applications, it is critical to have a controllable 
and efficient way to generate valley polarization of carri¬ 
ers. Although recent experiments have succeeded in gen¬ 
erating and modulating the valley polarization dynam¬ 
ically by circularly-polarized lights in monolayer transi¬ 
tion metal dichalcogenides, 9-12 it is yet desired to have 
valley control by static means in transport. 13 15 To this 
end, a valley filter is an indispensable component in val¬ 
ley tronics devices. Several schemes of valley filters have 
been suggested in the context of graphene, e.g., by using 
nano-constriction, 1 line defects, 16 strain engineering, 17 or 
the recently proposed valley-polarized quantum anoma¬ 
lous Hall phase. 5,18,19 Since in these filters propagating 
modes of both valleys generally co-exist, it is arduous 
to achieve perfect filtering due to intervalley scattering 
caused by lattice-scale defects, which are challenging to 
completely eliminate during fabrication. 

Here we propose a general scheme for creating a perfect 
valley filter in honeycomb-lattice materials utilizing a do¬ 
main wall (DW) 20 28 between two topologically-distinct 
insulating regions: a quantum anomalous Hall (QAH) in¬ 
sulator 5,29,30 breaking time-reversal symmetry (T) and 
a quantum valley Hall (QVH) insulator 3,5 breaking in¬ 
version symmetry (fP). Since both symmetries are de¬ 
stroyed, the valley filtering becomes efficient. It turns out 
that subgap-propagating chiral channels appear in the 
DW, which are dictated by the change of bulk topologi¬ 
cal charge across the DW. 20,24,31 Remarkably enough, all 
the unidirectional channels reside in the same valley when 
the valley-projected topological charge only changes at 
one designed valley across the DW. As a result, carri¬ 
ers in the filter could achieve a 100% valley polarization 
without any backscattering. Due to its topological origin, 
the performance of such a filter is robust against scatter¬ 
ing and structural imperfections. We show that both the 
valley index and the chirality of our proposed elementary 


filter can be easily controlled by external fields and can 
be probed by optical circular dichroism, which allows ver¬ 
satile valleytronics functionalities achieved in a relatively 
simple design. 

Toy model .— We first present our basic idea through 
the analysis of a simple toy model. Consider a spinless 
fermion model defined on a 2D honeycomb lattice (see 
Fig. 1(a)) 
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Here the first term represents the nearest neighbor hop¬ 
ping. The second term, first introduced by Haldane 29 
connects the next-nearest neighbors, with = + (—) 
if the electron makes a left (right) turn when hopping 
from site j to site i. 32 The last term is a staggered po¬ 
tential with opposite signs (£* = ±1) on different sublat¬ 
tices. 32 In the absence of the last two terms in (1), it is 
known that near half filling the energy dispersion is gap¬ 
less and linear down to the K and K' Dirac points, i.e., 
the two inequivalent corners of the hexagonal Brillouin 
zone. The two valleys are related by T and V symme¬ 
try operations. The gapless Dirac point at each valley 
is protected by a sublattice symmetry and the winding 
number of an infinitesimal Fermi circle. Both the Hal¬ 
dane term and the staggered potential open energy gaps 
at the Dirac points, as each breaks the sublattice sym¬ 
metry. Importantly, the Haldane term breaks T and the 
staggered potential breaks V. We consider a DW formed 
by a spatial separation of the two topologically-distinct 
states, e.g., 1 2 (A v ) to be nonzero only for region y > 0 
(y < 0), as shown in Fig. 1(b). Figure 1(c) plots the en¬ 
ergy spectrum with such a DW along the zigzag direction 
in order for the two valleys to be distinguished. In the 
bulk band gap, remarkably, there appears a chiral gap¬ 
less channel attached to valley K , localized at the DW, 
and propagating along +x direction. 

To gain a better understanding of this chiral gap¬ 
less valley-filtered channel, we focus on the low-energy 
physics of the toy model. Away from the DW, tak¬ 
ing t as the largest energy scale, we expand the Hamil¬ 
tonian around K (K') point and obtain a continuum 
model iL e ff = Hq + H/\, where the band term Hq = 


2 




y/a 


FIG. 1. The realization of a perfect valley filter, (a) The 
lattice model described by Eq. (1). Blue and yellow dots de¬ 
note the two sublattices. Representative next-nearest neigh¬ 
bor hopping are indicated by the red (green) arrows with 
h>ij ist + (—). (b) Schematics of the topological DW proposed 
in the main text, (c) The energy spectrum 33 of (b) with a 
valley-filtered chiral DW channel, (d) The valley-projected 
topological charges on the two sides of the DW. (e) Wave- 
function distribution for the zero energy mode at valley K. 
Here we take two different DW widths L (in units of a ) for 
a DW with a tanh (y/L) profile for the mass term. The solid 
lines are from Eq.(2) and the circles are solutions from the 
tight-binding model. 


VF{r z k x (Jx + kyCFy) and the mass term Ha is —A a T z cr z 
(X a = SVSt 2 ) for y -T Too and X v cr z for y -T — oc. Here 
the wave vector k is measured from the K ( K ') point, 
vp = 3at/2 is the Fermi velocity with a the distance be¬ 
tween nearest neighbors, r z = db denotes the two valleys, 
Pauli matrices a act on the sublattice pseudospin, and 
we assume that X v ,X a ,VF > 0. One observes that H e ff of 
valley K is a gapped two-band model with a mass sign 
reversal across the DW. It is quintessential that such a 
Jackiw-Rebbi problem 34 bears a zero mode: 
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where A is a normalization constant and the mass A 
is defined through Ha = A (y)cr z . This mode is expo¬ 
nentially localized at the DW, with a chiral dispersion 
E = vpk x around zero energy. In Fig. 1(e), we plot the 


wavefunction distribution of the DW zero energy mode 
for a sharp DW and for a DW with finite width. Benefi¬ 
cially, the presence of such a chiral gapless DW channel 
depends only on the designed signs of A at y = Too, re¬ 
gardless of the detailed DW characteristics, which reflects 
its topological nature. 

The robustness of the chiral gapless DW channel is 
intimately connected to the change of bulk topologi¬ 
cal charges across the DW. 5,20,24,31 Either side of the 
DW has a valley-projected topological charge J\f = 
±fd 2 kf2(k ); 20,22 where ft(fc) = V fc x A{k) and 
J4(fe) = {uk\iS7k\uk) are the Berry curvature and the 
Berry connection of the valence bands with \uk) being 
the periodic part of the Bloch eigenstate. The Berry cur¬ 
vature is out-of-plane for the 2D case and concentrated 
around each valley center in our model. For the con¬ 
tinuum model we obtain M = r 2 sgn(A)/2. Since 
A = —A a r z in the QAH domain (y > 0), J\f = —1/2 for 
both valleys, where the valley independence is dictated 
by V. Evidently, the total topological charge of both val¬ 
leys (Chern number) is one, featuring a QAH state. In 
the QVH domain with A = X v (y < 0), J\f = r^/2, where 
the opposition of topological charges of different valleys is 
a consequence of T (in this sense the QVH state is analo¬ 
gous to the quantum spin Hall state 5,32 ). These topologi¬ 
cal charges are illustrated in Fig. 1(d). The number of the 
valley-dependent gapless chiral DW modes is related to 
the difference of topological charges across the DW by an 
index theorem. 20,21,31 More specifically, for each valley, z/, 
the number of chiral modes moving in +x direction mi¬ 
nus the number of chiral modes moving in — x direction, 
is equal to the difference of the bulk topological charges 
J\f on the two sides of the DW, v = N(y < 0) —N(y > 0). 
It follows that v = 1 for valley K and v = 0 for valley K' 
in our case. As a consequence, one chiral gapless chan¬ 
nel is protected at valley K whereas none is at valley K '. 
This topological argument is consistent with our pervious 
analysis. 

The above analysis demonstrates that when across the 
DW the valley-projected topological charge changes (by 
an integer) for only one valley, the DW is guaranteed to 
possess valley-filtered chiral modes. Because all the prop¬ 
agating channels in the bulk gap are of the same valley, 
scattering cannot change the carriers’ valley index; and 
because of the chiral nature of these channels, there is 
no backscattering near the Fermi energy (assumed in the 
bulk gap). As a result, such a DW is perfectly transpar¬ 
ent for one valley but completely opaque for the other val¬ 
ley, hence behaving as a perfect valley filter. Regarding 
the valley-filter operation, the DW is hardly necessary to 
be exactly along the zigzag direction, as long as different 
valleys can be distinguished. Based on the topological 
argument, these DW channels are insensitive to the de¬ 
tailed DW profile, e.g., the orientation, and as long as 
the DW is smooth at the lattice scale and the intervalley 
scattering way below the Fermi surface is also suppressed. 
As one can observe in Fig. 1(e), the DW mode typically 
spans more than few tens of lattice sites even for a nar- 
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row wall and hence any scattering involving transfer of a 
sizable momentum becomes virtually impossible. 

Field control of valley filter. — Now we examine a more 
realistic system where spin is included and demonstrate 
the versatile operating modes of the valley filter which 
can be controlled by external electric and magnetic fields. 
Consider a tight-binding model Ho of a honeycomb lat¬ 
tice with the usual nearest neighbor hopping, as the t 
term in Eq. (1). Again, the QVH domain can be realized 
by a H-breaking staggered sublattice potential H v , as the 
X v term in Eq. (1). Physically, such a staggered potential 
can arise from the interaction with a substrate 35 or from 
the occupation of the two sublattices by different atoms 
or atomic orbitals. 7,36 Of particular interest is the case in 
which the two sublattices are slightly shifted with respect 
to each other out-of-plane forming a buckled structure. 
Such buckling has been found in 2D materials such as 
silicene, 37 germanene, 37 and quite a few chemically func¬ 
tionalized materials. 38 41 In this class of materials H v can 
be induced and tuned by a perpendicular electric field. 

Regarding the QAH domain, it can be engineered by a 
Rashba spin-orbit coupling H R = it R C L(*°^ x 

dij) z Cjp combined with an exchange coupling Hm = 
M'fZi a p c\ a s( fPcip, where dij is a unit vector pointing 
from site j to site i and s are the spin Pauli matrices. 
The T-breaking exchange coupling lifts the spin degener¬ 
acy producing a band inversion, whereas the ^-breaking 
Rashba term hybridizes the inverted bands yielding an 
energy gap. This QAH model has been proposed in the 
context of graphene. 30 The Rashba and the exchange cou¬ 
plings may be produced simultaneously by ferromagnetic 
adsorbed atoms or in proximity to a (anti)ferromagnetic 
substrate. 30,42,43 

A straightforward calculation yields that the valley¬ 
contrasting topological charge is r^sgn(A v ) for the QVH 
state with Hqvh = Ho + H v and sgn(M) for the QAH 
state with Hqah = Ho FH r + Hm, respectively. There¬ 
fore, at a DW between the two states, the number of 
chiral gapless DW channels at valley r z must be 

v = r z sgn(A„) - sgn(M), (3) 

based on the bulk-boundary correspondence as we dis¬ 
cussed in the previous section. Again the important con¬ 
sequence is that the chiral gapless DW channels must 
appear and only appear in a single valley. As an exam¬ 
ple, when X v > 0 and M < 0, v = 2 for valley K and 
o = 0 for valley K' . This implies there are two DW 
channels propagating in -\-x direction at valley K and no 
state at valley K\ which is confirmed by our numerical 
calculation and shown in Fig. 2(a). These DW states can 
also be solved in a low-energy continuum model, similar 
to what we have done for the toy model. 

It follows from Eq. (3) that both the valley index of 
the DW channels and their propagating directions (chi¬ 
rality) are determined by the signs of X v and M, which 
may be manipulated by tuning external electric and mag¬ 
netic fields. In Fig. 2(b), we schematically show how the 
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FIG. 2. (a) Energy spectrum of the monolayer model featur¬ 
ing two valley-filtered chiral DW channels (in red color). 44 (b) 
{ok,ok') as a function of the exchange coupling M and the 
sublattice staggered potential X v in the monolayer model, (c) 
(ok, vk') in the bilayer model, tuned by the interlayer poten¬ 
tial differences Uv and Ha on the two sides of the DW. 44 (d) 
Optical circular dichroism (with light frequency of the bulk 
gap) exhibiting contrasting feature across the DW. 


characters of DW channels vary as the signs of A v and M 
are switched. The two axes (X v = 0 and M = 0) corre¬ 
spond to topological phase boundaries at which the bulk 
gap vanishes. Eq. (3) and Fig. 2(b) are suggestive of a 
versatile control of the operating modes for the valley fil¬ 
ter. As an example, consider an initial mode with X v > 0 
and M < 0. Flipping only the sign of A v , the valley index 
becomes opposite whereas the chirality remains the same. 
Inverting both the signs of X v and M, the valley index 
maintains whereas the chirality is reversed. Both valley 
and chirality modes can be switched simultaneously by a 
sign reversal of M. All these characters can also be easily 
argued from the symmetry point of view. 

All the above discussions can be directly extended 
to a bilayer system that may exhibit a richer topologi¬ 
cal phase diagram due to the additional layer degree of 
freedom. We consider an AB-stacked bilayer graphene 
model following Ref. 45. Each layer is described by the 
same Hamiltonian Ho + H R + Hm as discussed above. 
The two layers are coupled by the interlayer hopping, 
Kt = t±J2(ie(t,A),je(b,B))4 c j > between the A site of 
top layer and the B site of bottom layer that are the 
nearest interlayer neighbors. Through dual gating, it is 
convenient to introduce an interlayer potential difference 
Hu = U JA Ci c i°i with Q = =b denoting the top and bot¬ 
tom layers. Hu plays a similar role here to that of the 
staggered sublattice potential H v in the buckled mono- 
layer system. Depending on the competition between the 
exchange coupling Hm and other potentials, various dis¬ 
tinct QVH and QAH states can be realized. 19,45,46 
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Compared with the monolayer case, of particular inter¬ 
est in the bilayer case is that its valley-projected topo¬ 
logical charge can have a magnitude of either one or two, 
depending on the model parameters. 45 As a consequence, 
at the DW between the bilayer QAH and QVH states, 
there exist two types of perfect valley filters with dif¬ 
ferent number of chiral DW channels: for one type it 
has two perfect valley-filtered channels, whereas for the 
other type it has four. In Fig. 2(c), we plot (z vk'\ the 
number of chiral channels in each valley, as a function of 
U a and Uy with other parameters fixed. Here we have 
assumed that the interlayer potential differences on the 
two sides of the DW, denoted by U a {Uy) for the QAH 
(QVH) domain, can be controlled independently. One 
notes that in the chosen parameter range, there are four 
regions of (z VK f ), of which (2, 0) and (4, 0) are perfect 
valley filters whereas the other two have DW channels in 
both valleys. Nonetheless, this indicates that by control¬ 
ling the interlayer potentials, it is possible to tune the 
conductance of a perfect valley filter, e.g., from 2 e 2 /h to 
4e 2 //z when the phase switches from (2,0) to (4,0). 

Optical circular dichroism . — Conventionally, gapless 
DW modes can be probed by transport or scanning 
tunneling spectroscopy. Since a perfect valley filter re¬ 
quires the breaking of T, optical circular dichroism may 
also be employed to detect the presence, the location, 
and the characters of the topological DW bearing a 
perfect valley filter. We sketch the circular dichroism 
in Fig. 2(d) for our model and illustrate the physics 
below using the standard argument of Fermi’s golden 
rule. The coupling strength of an interband optical 
transition with <j± circularly-polarized light is given by 
M±{k) = M x ±iM y , where = (u%\W k ^H\u v k ) is 
the interband matrix element. For a valley-projected 
chiral two-band system 5 with a gap 2A, e.g., gapped 
monolayer or bilayer graphene, the absorbance 4,47 of 
the <j± light with frequency uj has a universal form 
oc (1 ± 2t z A/uj) 2 Q(uj — 2A). When the optical frequency 


is close to the band gap, at valley K ( K ') the inter¬ 
band transition occurs only for cr + (cr_) circular polar¬ 
ization. For the QVH state, the net circular dichroism 
rj = (|AJ+| 2 — |AJ-| 2 )/(|.A4+| 2 + \M-\ 2 ) vanishes, as re¬ 
quired by T. In contrast, for the QAH state breaking 
T, there is a strong circular dichroism. Therefore, such a 
difference between the two domains in optical absorbance 
or photoluminescence can be used to probe the presence 
and the location of the valley filter. Since \M±\ 2 also 
depend on the sign of A, the circular dichroism can also 
probe to the character of the valley filter that is sensitive 
to the sign of A. 

Discussion .— We have proposed a general scheme for 
creating perfect valley-filtered channels at the DW be¬ 
tween two topologically-distinct insulating regions. The 
ease of local field control of each filter will allow the con¬ 
struction of more complex valleytronic devices using mul¬ 
tiple DWs or a percolating DW network. For example, 
a valley valve can be constructed by combining two val¬ 
ley filters in series. Whether a current can pass through 
the valve depends on whether the two filters have the 
same valley index and the same chirality, which can be 
controlled by local electric and/or magnetic fields. Al¬ 
though it is not our focus here to investigate specific ma¬ 
terials for its realization, we do note that it is promising 
to obtain the desired topological DWs in 2D materials 
with hexagonal lattice symmetry such as single or mul¬ 
tilayer graphene, silicene, and chemically functionalized 
monolayer materials, through controlled doping, gating, 
and substrate effects. Our scheme can be generalized to 
materials with more than two valleys, 48 and a perfect 
valley filter can also be created whenever across a DW 
the valley-projected topological charge only changes at 
one designed valley. 
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